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ABSTRACT 

In conjunction with the test and evaluation 
of the position fixing capabilities of the 
Army Manpack Loran Receiver AN/PSN-6, an 
extensive series of time difference and 
signal amplitude measurements were made 
within a 100 km map grid square encompassing 
Fort Monmouth, New Jersey. The test loca- 
tion is within the coverage area of the East 
Coast Loran-C Chain. The data were used to 
develop a simple "smooth-earth” model for 
the test area as well as to estimate the 
magnitude and distributions of deviations 
from this model. Local propagation pro- 
cesses associated with topographic features 
and the grid of overhead wires in the test 
area are shown to contribute to the devia- 
tions from the model. 

INTRODUCTION 

As part of a broad program to develop a capability for navi- 
gation and position fixing, the Army is in the process of 
developing a manpack loran receiver. The position fixing 
function of the receiver is to provide a real time display 
of either the loran time difference coordinates or the geo- 
detic coordinates of the receiver position. This conversion 
between the time difference and geodetic coordinate systems 
is accomplished by a small computer within the receiver. 

The research described here was designed to provide a data 
base for development of simple conversion algorithms as well 
as to provide an error budget for the resulting conversion. 
Since a fundamental variable of a loran system is the prop- 
agation time of the 100 kHz signal and since the manpack 
loran receiver is designed to operate on the ground, the 
present research was planned to provide information on 
ground wave propagation at the surface of the earth. These 
propagation results are the subject of this paper. 
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BACKGROUND 

The fundamental equation describing the functional relation- 
ship between a loran time difference and the loran chain 
parameters is 


TD(P) =[(D S - DJ/C] + ED, (1) 

where TD(P) is the loran time difference at a field point P, 
D is the great circle distance from the slave transmitter 
to the point P, D m is the great circle distance from the 

master transmitter to the point P, C is the propagation 
velocity of the loran signal, and ED is the emission delay, 
that is the sum of the propagation time from the master to 
the slave transmitter and the coding delay introduced at 
the slave station. [See Footnote (a).] 

The non-constant propagation velocity, which varies with 
the density and amplitude of terrain features and the 
electrical properties of the overland path, severely limits 
the utility of Eq. (1) for time difference estimation. On 
the other hand, the ground wave propagation velocity over 
sea water is a well-known quantity, so Eq. (1) is very use- 
ful for this application. A complete knowledge of the 
propagation velocity for all propagation paths is necessary 
for rigorous use of Eq. (1). This procedure requires an 
extremely large volume of data and clearly is not practical 
for a manpack loran receiver. This fact furnished the im- 
petus for the development of simplified coordinate conver- 
sion algorithms. In essence, the approach was to develop 
a simple local conversion model based on calibration and 
to investigate the accuracy characteristics of that model. 
The reader is referred to the work of Johler [1] and to the 
references cited therein for details on the complete treat- 
ment of overland loran propagation. 

EXPERIMENTAL PROCEDURE 

The test area in New Jersey is the 100 km square, 18T WV 
of the Universal Transverse Mercator Map System. A spec- 
ially equipped four-wheel drive mobile unit was used during 
data acquisition. Present instrumentation includes 2 mili- 
tary loran receivers, a timing receiver system with rubid- 
ium standard, and ancillary items including printers, 
oscilloscopes, and power supplies. The military receivers 
used roof -mounted whip antennas whereas the timing receiver 
used a rotatable roof-mounted loop antenna. 
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The calibration procedure is relatively simple, namely to 
obtain time difference readings at sites of known geodetic 
control. To obtain estimates of loran receiver performance^ 
approximately 100 time difference measurements were made 
with each receiver at each site. These measured values 
were then averaged to provide a time difference for each 
site. To eliminate separate site surveys, easily identifi- 
able topographic locations, such as road intersections, 
were used for geodetic control. Coordinates of all sites 
were determined from 7.5 minute USGS topographic maps. 
Criteria for site selection were positive identification, 
freedom from strong electromagnetic scatterers, and access- 
ibility with the mobile unit. In general, the absence of 
power lines was the most difficult criterion to meet. Sev- 
eral of the sites were at geodetic bench marks which pro- 
vided a higher order geodetic control. 

The distribution of calibration sites is shown in Fig. 1. 
The numbers shown in this figure are primarily for site 
identification purposes, but are also related to the time 
of calibration. Data at sites identified with numbers less 
than 1000 were obtained in December 1972, whereas those 
with identifiers greater than 1000 were obtained in July 
1973. This latter study was designed to give an increased 
calibration density within a 60 km square located in the 
SE corner of the primary test area. In addition, the ab- 
sence of nearby power lines was of extreme importance. It 
is estimated that for these locations, there were no wires 
within one kilometer of the site. For the December 1972 
measurements there were no wires within 300 meters of a 
site. This consideration of the proximity of overhead 
wires led to the classification of first and second order 
TD data, as indicated in Fig. 1. 

In October 1973, the timing receiver was used to obtain 
field strength information at the sites with identifiers 
less than 1000. Time difference measurements were also ob- 
tained to check the repeatability of earlier measurements. 

A study of the absolute phase variations of the Loran-C 
transmissions was initiated in October 1974. The sites for 
a ray path to the SS7-Y, Nantucket transmitter are desig- 
nated as phase track points in Fig. 1. All measurements 
were made relative to the reference point 1480 in Fig. 1. 
This point is approximately 1 km from the coast. The other 
sites are at approximately 10 km intervals along the ray 
path. The experimental procedure was to initialize at the 
reference point, make phase measurements at other sites, 
and then close the traverse at the reference site. Each 
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series of measurements required from 6 to 8 hours. Closing 
the traverse at the reference site provided an estimate of 
the frequency offset of the rubidium standard. 

The test area is a segment of the coastal plain which is 
essentially devoid of pronounced terrain features and con- 
ductivity discontinuities. Thus no significant perturba- 
tions of the loran signals were expected from these sources. 

ANALYTICAL PROCEDURE 

The model adopted for data analysis is a modification of 
Eq. (1). The underlying assumption for this model is that 
for modest coverage areas, a constant overland propagation 
velocity will provide a useful approximation. The equation 
is : 


TD(P) = (ED + a) + B(D s - DJ + €(0 $ , 6J . (2) 

In this expression, a and 8 are arbitrary parameters to be 
determined by least squares analysis of the measured data. 
The correction function, is bo account for sea 

water paths at a bearing angle 8 from the slave and master 
transmitters to the field point. This function is included 
to account for the significant difference between ground 
wave propagation velocity over land and over sea water. 

The parameter a can be interpreted as an average time dif- 
ference offset characteristic of the test area. The parame- 
ter 6 can be interpreted as the reciprocal of the local 
propagation velocity. However, in view of the simplicity 
of the model and the statistical method of analysis, strict 
physical interpretation of these parameters should be ap- 
proached with caution. [See Footnote (b).] 

For data processing, the variables D g and D m were calculated 

using the method of Sodano and Robinson [2] for the Clarke 
1866 spheroid [3], 

The correction e(8 ,0 ) was constructed from tabulated 

s m 

functions of the sea water path length as a function of 
bearing angle from the Nantucket and Carolina Beach trans- 
mitters. The sea water path functions were prepared from 
maps in increments of 5° in the bearing angle. The path 
length for intermediate bearing angles was determined by 
linear interpolation. The sea water path length functions 
are shown in Figs. 2 and 3 for the Carolina Beach and 
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Nantucket transmissions. In view of the simple interpola- 
tion procedure, the error associated with path length de- 
termination will be greatest for bearing angles greater 
than 25° for the Carolina Beach transmission. For this 
reason, all points with Carolina Beach bearing angles 
greater than 25° were not considered in the data analysis. 
This selection process left 61 calibration sites within 
the 100 km square. 

The correction function is 

e(0.0j = k[L (6 ) - L (0 ) ] . (3) 

s m s s m m v 

In this expression, L (0 ) is the sea water path length at 

.W 3 

a bearing angle 0„ from the slave transmitter and L (0 ) is 

the sea water path length at a bearing angle 0 m from the 

master transmitter, k = [(1/cp - (1/C )] where C^ is the 

ground wave propagation velocity over land and C g is the 

ground wave propagation velocity over sea water. The 

values used for these constants are: (1/C = 3.3416 ys/km 

and (1/C g ) = 3.3384 ys/km. Thus the constant k has the 

value 0.0032 ys/km. The value of C g was obtained from 

the tables of Johler and Berry [4] . The value of C^ was 

estimated from calibration of a geologically similar area 
in North Carolina where no sea water correction was required 

ANALYTICAL RESULTS 

As discussed previously, the data processing was designed 
to give the parameters a and 6 for the Nantucket and Dana 
slave configurations by the method of least squares. To 
evaluate the effect of coverage area size, the data were 
treated in two sets. One set included the total of 61 
points. The second set included the points within a 60 km 
square in the SE corner of the 100 km square. 

The analytical results are shown in Tables I and II. Also 
shown are the RMS deviations of the least squares fit for 
each data set. 
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TABLE I. STATISTICAL PARAMETERS FOR 100 km SQUARE 

(61 Data Points) 


Transmission 

a 

8 

RMS Deviation 

Pair 

(ys) 

(ys/km) 

for Set 

SS7-Y 

0.85 

3.346 

0.33 

SS7-Z 

-0.19 

3.339 

0.33 

TABLE II. STATISTICAL 

PARAMETERS 

FOR 60 km SQUARE 



(30 Data Points) 

Transmission 

a 

6 

RMS Deviation 

Pair 

(ys) 

(ys/km) 

for Set 

SS7-Y 

0.74 

3.346 

0.22 

SS7-Z 

-4.11 

3.349 

0.25 


An estimate of the experimental uncertainties was obtained 
by statistical analysis of the data acquired at each site. 
This procedure yielded an average value of 0.15 ys attrib- 
utable to instrumental jitter. In addition, it has been 
estimated that the use, of topographic maps introduces a 
location uncertainty of the order of 20 meters. For the 
test area, this corresponds to a time difference error of 
about 0.1 ys. Therefore, experimental processes are esti- 
mated to contribute an uncertainty of the order of 0.18 ys. 

The distribution of the magnitude of the time difference 
deviations is shown in Figs. 4 and 5 for the Nantucket- 
Carolina Beach (SS7-Y) and the Dana-Carolina Beach (SS7-Z) 
configurations, respectively. The solid curves are the 
normal distributions corresponding' to the standard devia- 
tions calculated for each slave configuration. The areal 
distribution of the time difference deviations for each 
slave configuration is shown in Figs. 6 and 7. 

The results of the field strength study of October 1973 are 
shown in the contour plots of Figs. 8 and 9 for the Nan- 
tucket and Dana transmissions, respectively. The contour 
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Fig. 4. SS7-Y Time Difference Deviation Magnitude Distribution. 
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Fig. 5. SS7-Z Time Difference Deviation Magnitude Distribution. 
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Fig. 6. Map Showing SS7-Y Time Difference 
Deviations from Theoretical Model. 


328 







Fig. 7. Map Showing SS7-Z Time Difference 
Deviations from Theoretical Model. 
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values are the output voltage of the amplitude strobe of 
the timing receiver when tracking the third cycle crossover 
with an input attenuation of 40 dB. The indicated conver- 
sion constant of 46 yv/m/volt was estimated from the an- 
tenna characteristics to provide a corresponding approxi- 
mate value of the field strength. 

The phase of the Nantucket transmission relative to the 
reference point as a function of the distance from the 
transmitter is shown in Fig. 10. The plotted points are 
an average of measurements on six different days. The devi- 
ations of these measurements from a least squares fitted 
straight line are shown in Fig. 11. The experimental re- 
sults have been corrected for loran chain variations from 
data furnished by the United States Coast Guard. The error 
bars in Fig. 11 represent typical uncertainty estimates 
arising from rubidium standard frequency offset. 

DISCUSSION AND CONCLUSIONS 

The simplified mathematical model presented provides a 
reasonably accurate description for a small segment of the 
coverage area, namely the 100 km square area. The observed 
magnitude of the deviations from the model are randomly 
distributed. Specific parameters of the model are sensitive 
to the size of the coverage area. In view of the simplicity 
of the model, strict interpretation of the parameters in 
terms of propagation properties is not possible. For exam- 
ple, the value of the parameter B for the SS7-Z slave shown 
in Table I (3.339 ys/km) is essentially the value expected 
for an all sea water path. [See Footnote (c) . ] 

The time difference deviations and field strengths exhibit 
a pronounced areal variation. Furthermore ■, the contours 
of both variables exhibit a preferred orientation in a 
NE-SW direction. Figure 12 shows salient topographic and 
geological features of the test area. These features also 
exhibit a preferred orientation in a NE-SW direction. Con- 
sequently, there appears to be some correlation between 
contour orientation and the topographic and geological 
features . 

The standard deviations of the time difference data are 
0.33 and 0.22 ys for the 100 km square and the 60 km square, 
respectively. These values exceed the 0.18 ys estimated to 
arise from experimental sources. Clearly the results for 
the smaller test area are in better agreement with the 
theoretical model than are the results for the larger test 
area. The criterion for proximity of nearby overhead wires 
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also influences the deviations from the theoretical model. 
Although the time difference data does not allow a defini- 
tive separation of the contributions of area size and prox- 
imity of overhead wires, both factors, as well as previ- 
ously discussed topographic effects, appear to contribute 
to the deviations from the idealized model. 



Fig. 12. Map of New Jersey Showing Salient 
Topographic and Geologic Features. 


Preliminary measurements of absolute phase along a ray path 
from the Nantucket transmitter across the 100 km test area, 
yielded a linear variation of phase with distance from the 
transmitter. The least squares slope is 3.352 ys/km, with 
a standard deviation of 0.3 ys. These observations are 
considered to be consistent with the time difference meas- 
urements for the 100 km square. Since this particular ray 
path is located in a region of nearly constant conductivity 
and is devoid of terrain irregularities, it is assumed that 
the major contribution to the observed deviations arises 
from scattering associated with overhead wires. 
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These results are of importance to the position fixing 
accuracy of ground -deployed loran receivers. Randomly dis- 
tributed time difference deviations of the order of 0.3 us 
from an idealized model have been observed. This study 
suggests that scattering associated with topographic fea- 
tures as well as from man-made sources such as overhead 
wires contribute to the deviations. The results yield a 
realistic standard deviation of position location accuracy 
of the order of 60 meters. 
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Additional Background Informati on 

(a) The assumption of identical propagation velocities 
over the paths D g and D m is implicit in Eq. (1). 

(b) The interpretation of 3 as the reciprocal of a local 
propagation velocity is valid only to the extent that 
the propagation velocities over the paths D and D 
are identical. 

(c) Furthermore, a change in the parameter a of 4 ys for 
two fits of the same area is devoid of physical sig- 
nificance. This behavior is related to the fact that 
D g and D m are large quantities so that in the least 

squares process, small variations in the parameter 3 
are compensated for, by large variations in the parame 
ter a. Additional physical constraints as discussed 
by Doherty [5] are necessary for a realistic physical 
interpretation of a statistical model. 
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